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ABSTRACT 

Aims. Stars in the narrow mass range of about 2.5 and 3.5Mo can develop a thermally unstable He-burning shell during its ignition 
phase. We study, from the point of view secular stability theory, these so called thermal micropulses and we investigate their properties; 
the thermal pulses constitute a convenient conceptual laboratory to look thoroughly into the physical properties of a helium-burning 
shell during the whole thermally pulsing episode. 

Methods. Linear stability analyses were performed on a large number of 3Mo star models at around the end of their core helium- 
burning and the beginning of the double-shell burning phase. The stellar models were not assumed to be in thermal equilibrium. 
Results. The thermal mircopulses, and we conjecture all other thermal pulse episodes encountered by shell-burning stars, can be 
understood as the nonlinear finite-amplitude realization of an oscillatory secular instability that prevails during the whole thermal 
pulsing episode. Hence, the cyclic nature of the thermal pulses can be traced back to a linear instability concept. 

Key words. Stars: AGB and post-AGB - Stars: interiors - Stars: oscillations 



1. Introduction 



Accor ding to the published literature, iGurevich & Lebedinskil 
(Il947h seem to have been the first to study the thermal sta- 
bility of nuclear-burning shells in stars; the authors aimed at re- 
vealing the physical nature of novae and supernovae. The work 
of Gurevich & Lebedinksi apparently went unnoticed outside 
of the Soviet Union so that 15 years passed before the field 
advanced: The seminal paper of iSchwarzschild & HarmI (Il965h 
reported the thermal instability of the He-burning shell in a 
IM0 asymptotic giant-branch (AGB) star The instability was 
stumbled upon once the time-step of the evolution computation 
was chosen sufficiently small during the pertine nt evolutionary 
phas e. Independently of Schwarzschild & Harm. [WeigertI d 19651 
11966 ) observed the same phenomenon at the beginning of the 
double-shell burning phase of a 5M0 star model. Weigert's com- 
putations in dicated already the cyclic nature of the instability. 
iRosd d 19661) . studying the advanced evolution of a O.53M0 he- 
lium star, also picked up thermal pulses once the evolutionary 



' The paper referred to in Zel'dovich & Novikov's "Stars and 
Relativity" is a "translation" for astronomers that appeared in 1955 
in Trudy chetvertogo sovshchaniya po voprosam kosmogonii, 143. The 
early work by Gurevich & Lebedinski seems to have made its way to the 
west essentially through the reference in the Zel'dovich Sl Novikov text, 
as whenever astronomers refer to Gurevich & Lebedinksi, they usually 
cite the 1955 version and not the technical original of 1947. 

^ The 5Me star studied by Weigert was previou sly computed into th e 
advanced evolutionary stage by Kippenhahn fsee'Kippe nhahnl l ll965h ). 
He could not get the star to eventually turn into a white dwarf as a 
carbon flash developed on the AGB. The star's center heated up too 
much as neutrino cooling was neglected. Weigert's computations in- 
cluded neutrino energy losses; this energy sink kept the central tem- 
perature low enough to prevent the onset of central carbon burning and 
the star evolved further along the AGB to finally run into the thermal 
instability of the He burning shell. 



time-step became sufficiently short. As Rose was interested in 
the evolution to the white-dwarf domain of He-star models, he 
chose large time-steps across the helium shell burning phase 
so that he captured only the last t hree thermal pulse cycles of 
the instability episode. In any case. iRos^ (11966) emphasized, as 
Weigert did, the cyclic nature of the phenomenon and he was 
first to present eigenfunctions in the suiToundings of the unsta- 
ble nuclear burning shell. 

Alreadv [Ledouxl (1 19621) noted that, in principle at least, the 
secular stability problem admits of complex eigensolutions. The 
st ability analysis of thin nuclear -burning shells was extended 
in iHarm & Schwarzschildl (Il972l) wherein the thermal pulses, 
eventually considered as a cyclic phenomenon, were associated 
with the star's secular modes. The authors connected the onset of 
the thermal pulses with an overstable oscillatory secular mode. 
After the first few pulse cycles, however,"[. . . ] the oscillatory in- 
stability gives way [. . . ] to a simple exponential instability which 
causes the repetitive helium-shell flashes described in earlier in- 
vestigations [. . . ]" as Harm & Schwarzschild ( 1972) put it. This 
last statement seems to have been taken as the final word on the 
issue and was perpetuated as such in the astronomical literature 
thereafter 

Despite the apparent settlement of the case, the point of view 
adopted bv lHarm & Schwarzschildl (1 1 9721) continued to intrigue 
as it does not explain the cause of the cyclic nature of the ther- 
mal pulses once the oscillatory instability gives way to an ex- 
ponential one. The numerical studies showed that the period of 
the initial oscillatory secular mode was close to the length of 
the developing thermal pulse cycles of the helium shell. Why or 
how the star remembered t he linear result so closely remained 
unexplained ( Hansen|[l978l) . If the cyclic nature of the thermal 
pulses is indeed a nonlinear phenomenon that escapes accessi- 
bility of a linear analysis, then the "relaxation oscillations" of 
the thermal-pulse phase are eventually one case of the realiza- 
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tion of hard self-excited o scillations (ICoxl 19801) in stars (see also 
iBuchler & Perdandll979l) . 

Besides the thermal pulses that are associated with secularly 
unstable thin helium-burning shells in low- and intermediate- 
mass stars along the asymptotic giant branch, secularly un- 
stable shell burning is also encountered at the onset of off- 
center helium burning in low-mass stars as they evolve from 
the tip of the red-giant branch onto the horizontal branch (e.g. 

oma sl ll967l:lDespaii3ll981h . But also more exotic episodes 
of stars' lives can get entangled in He-shell instabilities: X- 
ray bursts as observed in mass-acc r eting neutron stars (e.g. 
[Hansen & Van HornI 119751: 'Bildsten' Il995h . Even hydrogen- 
burning shells on accreting white dwarfs can go cyclical ly unsta- 
ble (e.g. lGiannone & Weigertlll967l:ICassisi et al.lll998l and ref- 
erences in the latter) via the same physical mechanism. Hence, 
a solid understanding of the thermal pulses is beneficial for dis- 
parate fields of stellar astrophysics. 

This paper deals with those thermal pulses that are en- 
countered when, under suitable circumstances, a helium shell 
builds up an d ignites at the end of co re He burning of a 
3M0 star (Mazzi telh & D'Antonal Il986l) . Applying the sim- 
ple, local instability cri t erion f o r thin nuclear-burning shells of 
ISchwarzschild & HarmI (Il965h . iMazzitelli & D'Antonal (Il986l) 
concluded that the instability was physical and due to the same 
mechanism as encountered later, higher up along the AGB. The 
thermal pulses at the onset of He shell burning were termed 
thermal micropulses (ThMPs) as the surface luminosity varies 
only marginally as compared to the variation during a thermal 
pulse on the advanced AGB. In an attempt to evolve a star 
from the main sequence to the terminal Debye-cooling stage as 
a white dwarf, therm al micropulses were also encountered by 
lAlthaus et al.l (|2002|) in stellar models for which semiconvection 
was neglected in the computations. It seems that only stars in the 
narrow mass range of 2.5 -3.5Mo develop ThMPs under suitable 
conditions. 

The evolutionary tracks of stars undergoin g ThMPs (cf. 
IMazzitelli & D'Antonall 1 9861: lAlthaus et al.ll200 2') show that the 
manifestation of the pulses in the stars' observables are mini- 
mal. The pulse amplitudes in the total luminosity are small and 
the time-scale of the pulses so short that eff'ects on a population 
of stars on the early second ascent of the giant branch are very 
unlikely to be detectable even in large stellar aggregates. Despite 
ThMPs being unlikely to be observable or possibly even realized 
in nature, their very existence in star models constitutes a wel- 
come opportunity: With the ThMPs, we can theoretically study 
the stability properties of a secularly unstable nuclear-burning 
shell from onset of the instability to the very end when the shell 
stabilizes itself again. In particular, the self-stabilization of the 
shell distinguishes the ThMPs from thermal pulses along the up- 
per AGB. Furthermore, we expect that the properties of the sec- 
ular mode spectrum of the ThMP phase are also representative 
for the thermally pulsing stars on the advanced AGB. 

In the following, we use the ThMP phase as a idealized lab- 
oratory to study the episode of cyclic thermal pulses once again 
from the point of view of linear secular stability theory. We take 
advantage of the result of Gabriel (1972) that over time-scales 
that are short compared with a star's nuclear time-scale its quasi- 
static evolution can be expanded in terms of secular eigendata. 

We follow the pertinent part of the secular eigenspectrum 
throughout the whole ThMP episode and we show that ThMPs 
can be understood as an oscillatory secular instability; hence 
nonlinear theory via relaxation oscillations does not have to 
be invoked to explain the existence of the cyclic instability. 
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Fig. 1. Evolutionary track of a 3M0 star with pertinent zoom-ins 
to the epochs of core He burning (panel b) and to the develop- 
ment of double-shell nuclear burning (panel c). Core He burning 
starts at the tip of the first giant branch that is seen in the up- 
per right of panel a. The transition from He burning in the core 
to He burning in a shell leads to a temporary luminosity drop 
(after model 1490 which is designated in panel c); during this 
luminosity decrease, the thermal micropulses set in. After the 
He-burning shell is well developed, the micropulses die out and 
the star continues its second ascend of the giant branch towards 
the asymptotic giant branch. The epochs of the begin and the end 
of the linear secular stability analyses are indicated by the letters 
A and E, respectively. The location of model numbers 1400 and 
1490 along the evolutionary track are given for later reference in 
the text. 

Nonlinear effects, on the other hand, determine the temporal 
shape of the observables of the pulses and they seem to govern 
the amplitude evolution of the pulse cycles. 



2. Modeling Procedure 

Stellar evol ution was computed with the LPCODE c ode de - 
scribed in lAlthaus et all (|2002|) and lAlthaus et all ( l2003h . 
Convective overshooting was modeled as an exponentially de- 
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caying diffusive process above and below every convective 
shell, including the boundary of the convective core (during the 
main-sequence and central heliu m-burning phase) an d the con- 
vective envelope (for details see lAlthaus et al.ll2003b . The nu- 
clear network considered for the stellar modeling accounts ex- 
plicitly for 16 chemical elements and 34 thermonuclear reac- 
tion rates to follow hydrogen and helium burning. Abundance 
changes of the nuclear species were included by means of a 
time-dependent scheme that simultaneously treats nuclear evo- 
lution and mixing processes. The efficiency of convective mix- 
ing was described by appropriate diffusion coefficients which 
are specified by the treatment of convection. The equation of 
state included partial ionization, ionic contributions, partially de- 
generate electrons, and Coulomb interactions. Radiative opac- 
ities, including carbon- an d oxygen-rich compositi ons, were 
those from the OPAL sets jlglesias & Roperslll996l) . comple- 
mented, at low temperatures, with the molecular opacities from 
[Alexander & Fergusoiil (Il994l) . 

The computations to investigate the linear secular stabil- 
ity properties of the star models were performed with the 
Riccati code for radial nonadiabatic stellar pulsations; the cur- 
rent version of the code is based on the one described in 
iGautschv & Glatzel (Il990l) . In addition to the canonical radial 
pulsation equations, we used also modified stability equations 
wherein the condition of thermal equilibrium of the background 
model was relaxed. For completeness, the relevant linear stabil- 
ity equations are presented in Appendix A. 



3. Results 

This section contains the results obtained from the computation 
of the evolution of a 3Mq model with initially homogeneously 
distributed abundances of X = 0.705, Y = 0.275. We followed 
the star's evolution from the main sequence to stable double- 
shell burning on the lower AGB. The second part of this section 
is devoted to the secular stability analyses performed on the evo- 
lutionary star models between about the ignition of the helium- 
burning shell and the termination of the thermal pulses. 



3.1. The micropulses as seen in stellar evolution 
computations 

Figure [T| shows, on three zoom levels, the evolution on the 
Hertzsprung-Russell (HR) plane of the 3Mq star mentioned be- 
fore. The displayed track resulted from the computations without 
convective overshooting - thermal micropulses developed only 
in such stellar models; this aspect will be returned to further 
down. Panel (a) captures the star's locus between core hydrogen 
burning on the main sequence and the establishing of two burn- 
ing shells at the bottom of the AGB. Core helium burning sets in 
at the tip of the first giant branch at the upper right of panel (a). 
Panel (b) shows the stump of a blue loop that develops during 
core helium burning. After passing through a local luminosity 
minimum at logL/L© = 1.792, with a central helium abundance 
Fc = 0.8503, the star's luminosity rises again to pass through 
a second local luminosity maximum (logL/L© = 2.015). Label 
A denotes the epoch of the evolution at which the secular sta- 
bility analyses were started. For ease of presentation, the sec- 
ular stability results will be displayed mostly as a function of 
model number; label A corresponds to model number 1300 with 
Fc = 0.288. The heHum shell ignited at = 0.04; the corre- 
sponding model number, 1400, is labeled along the track in panel 
(b). At the second local luminosity maximum (model number 
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Fig. 2. Variation of various luminosities during the thermal mi- 
cropulse phase of the evolution of a 3M0 model. Ltot stands for 
the total stellar luminosity; it is composed of the helium-burning 
one (Liie), the luminosity generated in the hydrogen shell 
(LcNo), and the luminosity induced by expansion/contraction 
(Lg) in the own gravitational field. The abscissa measures the 
age of the star in 10^ years (My). 



1490 in panel c), Y^ dropped to 0.0325 and the triple-o' energy- 
generation rate of the nuclear-burning shell reached the same 
strength as that of core helium burning. Hence, after model 1490, 
the star's helium burning is shell dominated. Finally, panel (c) 
zooms in to the cyclically varying luminosity and effective tem- 
perature that develop shortly after model 1490. The star passes 
through a local luminosity minimum after the 6th pulse cycle. 
Around this local luminosity minimum, the surface-luminosity 
variation reaches its maximum with a relative amplitude of 5 %. 
During the ensuing upward evolution along the lower AGB, 
the amplitude of the ThMPs decreases continuously (see also 
Fig.|2]i. The variability disappears completely at about the epoch 
labeled by E (model number 3234) when the helium in the center 
has essentially burned away, i.e. having reached Y^ = 1.15- lO^^*. 

The temporal evolution of pertinent luminosity components 
is shown in Fig. |2l At model 1490, the epoch of the onset of the 
ThMPs, the star has reached an age of about 395 My. The am- 
plitude of the luminosity variation in the He-burning shell (Lhb) 
grows rapidly to reach maximum amplitude - varying between 
10 and 1 IOL0 - in cycle 7. Since the layers in and above the He 
shell expand during the increase of the shell's energy produc- 
tion, the hydrogen shell effectively cools so that its luminosity, 
LcNO, drops. Figure |2] shows also the mirror-image behavior of 
the mechanical work, Lg, in the star's gravitational field. Energy 
is required for expansion during the Lee rise and it is gained dur- 
ing the shrinking of the envelope around Lee minimum. In the 
sum of the partial luminosities, Ltot = ^He + ^cno + L^,, the total 
luminosity variation remained always very small, less than 5 %. 
It is mainly Lg which counteracts the strong variation of Lhb, de- 
spite the fact that hydrogen burning mostly dominates the energy 
production. 

Figure 3 shows, for the H-burning shell (top panel) and the 
He-burning shell (bottom panel), the temporal evolution of phys- 
ical quantities at the center of the corresponding nuclear-burning 
shell over the first nine ThMPs. The data for temperature, T, den- 
sity, p, and the fractional mass, q, were shifted vertically to fit 
into a decently scaled single diagram. The vertical dashed lines 
in Fig. 3 mark the epochs of maximum Lhe- The mean radii of 
the H- and He-shell slowly increase with time reflecting the slow 
ascent of the stellar model along the second giant branch (cf. 
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Fig. 3. Temporal variation of selected physical quantities at the 
center of the hydrogen-burning shell (top panel) and the helium- 
burning shell (bottom panel). The vertical dashed lines indicate 
the epochs of the first nine Lhb peaks of the ThMP episode. The 
loci of temperature (T), fractional mass (q), and density (p) were 
shifted vertically to fit the figures that were set up to display the 
variation of the radii (r) of the shell centers. 

Figure[TJ. Notice the quick saturation of the amplitude growth of 
the radius variation of the He shell after about the first four pulse 
cycles. The increase in q, the fractional mass, of both shells is 
small. In the He-shell, the density minimum occurs at minimum 
radius of the shell and vice versa; maximum temperature is well 
in phase with Lee reflecting the dominating dependence on tem- 
perature of the nuclear energy-generation rate of helium burning. 
The H-shell behaves differently: Minimum density agrees with 
minimum radius and with maximum temperature (this coincides 
with the maximum-density phase of the He-shell). Hence, mini- 
mum radii of H- and He-shell coincide whereas maximum radius 
of the H-shell is reached slightly later than maximum radius of 
the He-shell. The wiggling in log p at the center of the H-shell 
is an artifact contracted from our definition of the shell center 
as the gridpoint with maximum nuclear energy generation rate; 
this maximum slightly shifts in mass and as the density contrast 
across the H-shell is significant, numerical noise is easily picked 
up. 

After a total of 30 pulse cycles, at about 403.7 My (model 
number 3234), the ThMPs die out. Hence, the ThMP episode 
lasts for about 8.7 My. The mean pulse-cycle length (tthmp) is 
0.29 My, a value that remains very stable throughout the whole 
instability episode. 

3.2. The secular stability analyses 

Between model numbers 1300 and 3234, i.e. between labels 
A and E of Fig. [T] we performed linear stability computations 



on the evolutionary models to study their secular eigenspectra. 
The physical quantities were assumed to change proportional to 
exp(i<ijf). The quantity cr that is used in the following is the fre- 
quency o) divided by sfiGMJR^ with G being the gravitational 
constant. To obtain approximate e-folding times in mega-years 
in Fig.m compute 5.12 x lO^^/cri. 

For the stability computations we relaxed the requirement 
of thermal equilibrium of the stellar models; the assumption of 
complete equilibrium is what is usually implemented in stel- 
lar pulsation analyses. For the secular problem, in particular 
when pulses develop, Lg contributions might be significant so 
that refraining from the assumption of thermal equilibrium ap- 
pears more appropriate. The linear stability equations, contain- 
ing approximations for the thermal imbalance terms, are given 
in Appendix A. 

The left panel of Fig. |4]shows cri, the imaginary parts, of the 
lowest few secular eigenmodes as a function of model number. 
Complex conjugate eigensolutions (oscillatory secular modes) 
are indicated by full dots. Dot-free lines indicate vanishing real- 
parts of cr, i.e. purely monotonic modes. For better correlation of 
the linear results with the evolutionary computations, the right- 
hand panel of Fig.|4]plots the Lee oscillations. For completeness, 
the age in mega-years is labeled along the rightmost ordinate of 
the figure. 

The complexity of the secular eigenspectrum is consider- 
able; and this is so already before the onset of the ThMPs. 
During the late core helium-burning phase, the modal diagram 
is a complicated web of complex-conjugate modes unfolding 
into pairs of non-oscillatory modes that later merge again into 
complex-conjugate solutions in different arrangements. In all 
circumstances, the number of modes involved is conserved. 

At model 1300, two unstable monotonic secular modes were 
encountered. As the core dominance of the nuclear burning 
ceased, the strength of the instability diminished. Eventually, the 
most unstable secular modes merged with other branches into 
oscillatory modes. Most of them became or remained stable dur- 
ing the early helium shell-burning phase. A discussion of the 
secular-stability behavior during the core He-burning phase of 
intermediate-mass stars is postponed to a forthcoming paper 

The onset of the ThMP phase manifested itself in the least 
damped oscillatory secular mode going overstable. The first two 
ThMP cycles appeared as unstable oscillatory secular modes in 
the linear analysis. Beginning with cycle three, the unstable com- 
plex conjugate eigensolution unfolded into a pair of monotonic 
modes during the maximum of Lhb within a cycle. During phases 
of reduced efficiency of the helium burning shell, the two mono- 
tonic modes merged to go into the next cycle through a short 
damped episode. The cyclic unfolding of the unstable oscillatory 
secular mode continued to cycle 21 after which the amplitude of 
the ThMP was seemingly small enough for the unstable secular 
mode to remain purely oscillatory. The ThMPs, as seen in Lee 
(e.g. right panel of Fig. [Hi, terminate by the dominant unstable 
oscillatory secular mode going stable but remaining oscillatory. 

The details of the behavior of the complex eigenvalue cr 
throughout the phases of the ThMP episode must not be over- 
interpreted. After all, the stability analyses were performed on 
the thermally pulsing models themselves. Still, to gain access 



^ In analogy with the pulsation theory of stars, our linear secu- 
lar analyses can be compared to performing linear pulsation analyses 
on nonlinearly pulsating star models. The period of pulsation derived 
from such a linear analysis varies cyclically proportional to R,{tf^^ 
where Rt{t) is the nonlinearly pulsating stellar radius. Also the damp- 
ing/excitation rate (our a\ here) varies along a pulsation cycle. In par- 
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Fig. 4. The left panel of the graph shows the modal 
diagram with the imaginary parts (cri) of the secular 
modes of the model series presented in Fig.[T] The evo- 
lutionary epochs defining the ordinate are counted in 
model numbers on the left and in mega-years on the 
right-hand side. Oscillatory secular modes, with non- 
vanishing real parts of the eigenfrequencies, are dis- 
played as filled circles. For reference in the text, the 
three main branches of oscillatory secular modes ex- 
tending over most of the ThMP phase are labeled by 
Oi, O2, and O3. To ease the association of model num- 
bers with a characteristic physical quantity of that evo- 
lutionary epoch, the panel on the right displays the tem- 
poral evolution of the luminosity generated by helium 
burning. 



to the conditions under which complex secular eigenvalues un- 
fold into real pairs, secular analyses on the full cycles of ThMP 
models might prove helpful. 

In addition to the unstable, mainly oscillatory secular mode, 
Fig.Hlshows a damped oscillatory mode, O2 that never unfolded 
during the whole ThMP episode. The oscillatory secular mode 
baptized O3, on the other hand, unfolded again cyclically into 
pairs of monotonic modes, but in this case during the minimum 
phases of Lhb; Fig.|4]shows only one of the monotonic branches 



ticular when the nonlinearly pulsating model is hottest, the pulsational 
excitation by the %-mechanism of the stellar envelope is diminished as 
the driving zone lies too close to the surface; at the phase of lowest sur- 
face temperature, excitation is diminished as the driving region lies too 
deep in the envelope. 



of O3 since a follow-up of both of them turned out to be compu- 
tationally cumbersome due to other close-lying secular modes. 

In the range of the tri domain considered here, we en- 
countered only three oscillatory or partially oscillatory secular 
modes. At least in the cri - model number (or age) diagram these 
oscillatory modes overlap with purely monotonic secular modes. 
Except for the mode with the longest time-scale, all of the mono- 
tonic modes are always damped. Note that the crossing of modes 
in Fig.|4]is a projection effect due to plotting cri only instead of 
the complex cr; on the complex cr plane, the modes are well sep- 
arated. 

The secular mode with the smallest cri values and which re- 
mains always monotonic switches between weakly stable and 
unstable; the growth rates of this monotonic mode are two orders 
of magnitude smaller than those of the overstable complex secu- 
lar mode. After model 3216, when the oscillatory secular mode 
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went stable and ThMPs disappeared, the lowest monotonic sec- 
ular mode remained unstable with an e-folding time of 3.4 My 
(as a comparison, an evolutionary change of the total luminosity 
by 5% takes 3.1 My at this epoch of evolution). 

For the set of secular modes that remain monotonic at all 
epochs after the onset of the ThMPs, the ixi's vary all in phase 
during a pulse cycle. Furthermore, the general ctj evolution dur- 
ing the ThMP episode is the same for all these monotonic secu- 
lar modes; and this evolution differs from the oscillatory secular 
modes. The latter fact is most evident in the evolution of cri of 
the most strongly damped oscillatory mode relative to the mono- 
tonic secular modes between models 1500 and about 1800. 

The stellar model sequence whose computation included dif- 
fusive overshooting and which showed no evidence at all of ther- 
mal pulses in the evolution computations also revealed no over- 
stable oscillatory secular mode at any time during the double- 
shell nuclear burning phase. However, from the very beginning 
of the core helium burning through the double-shell burning 
phase at least one monotonic secular mode was always unstable. 
Overshooting left its imprints in the abundance structure and in 
a bigger stellar core as compared with ThMP models. 

3.3. On the suppression of thermal micropulses 

The ThMPs are found to be a delicate phenomenon. First of 
all, ThMPs are encountered only in computations with appro- 
priate convection treatment and in the narrow mass range be- 
tween about 2.5 and about 3.5Mo. The easiest and in some 
sense the most dangerous way to miss the ThMPs is by 
means of the time-step choice to evolve a model star. Already 
ISchwarzschild & HarmI (11965) cautioned in their paper that in- 
appropriately large time-steps, of the order of twice the e-folding 
time of the instability, suffice to numerically suppress a pulse. 
We repeated the computation described in Sect. 13.11 (we will 
refer to it to as the "original" one) by enforcing a larger time- 
step. The coarser-grained time evolution sequence started helium 
shell burning with a time-step that was 1 .88 times larger than the 
one of the original sequence. The structure (position and thick- 
ness) of the helium shell and the abundance profiles differed only 
marginally from those in the original computation. Nonetheless, 
the helium shell remained essentially stable: We only observed 
a short initial oscillatory phase of Lhb lasting for roughly three 
strongly damped cycles. The first cycle had the largest ampli- 
tude with a relative Lhb variation of 11 %. (Remember, in the 
original evolution computation, the relative ThMP amplitude of 
the helium shell reached 400 %.) No variability was manifested 
in Ltot- Even an inspection at highest numerical resolution of the 
computed observable stellar parameters revealed no sign of an 
instability whatsoever (see panel c of Fig.|5]). 

Figure |5] displays results from computations performed on 
the pulse-suppressing stellar evolution sequence. As mentioned 
before, we found one overstable secular mode with characteris- 
tics that were comparable to the overstable secular mode in the 
thermally pulsing model sequence. The dashed horizontal lines 
extending across all three subpanels of the figure indicate, at the 
bottom, the onset of the secular instability and the termination 
of the instability in the upper part of the plot. Panel (a) shows 
the period evolution of the overstable secular mode as a func- 
tion of model age; the dash-dotted line in this panel indicates the 
mean pulse period as derived from the fully pulsing evolution- 
ary sequence. Note that the period obtained in the linear analysis 
deviates by less than 10% from the completely independently 
computed nonlinear cycle length of the pulses. The inverse of 
the e-folding time (measured in mega-years) of the instability 
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Fig. 5. Panels (a) and (b) show eigendata of the overstable oscil- 
latory secular mode as computed in the pulse-suppressing model 
sequence which is characterized in panel (c) by the helium- 
luminosity (Lue) and the total luminosity (Ltot) evolution. Panel 
(a) displays the period evolution of the overstable oscillatory sec- 
ular mode as a function of model age. The vertical line shows the 
magnitude of the average pulse-cycle length that was deduced 
from Fig. |2] The inverse of the e-folding time Te (measured in 
units of mega-years) of the overstable secular mode is plotted in 
panel (b). 

is plotted in panel (b). For most of the instability phase, the e- 
folding time, re, is less than or of the order of the evolutionary 
time-scale itself. Panel (c) finally shows the luminosity evolu- 
tion of helium burning (Lhc) and of the star's total luminosity, 
Ltot- The quantity Lhc of the evolution computation shows only 
a weak sign of instability that lasts for roughly three strongly 
damped cycles; the total luminosity carries even no sign at all 
of an instability. We emphasize again that the phase of instabil- 
ity as obtained from the linear secular analysis is comparable to 
the evolutionary phase over which the ThMPs are encountered 
in the "non-suppressing" model sequence; this is traced out with 
the vertical line plotted from the local maximum of Ltot in panel 
(c). The local maximum at the bottom of the figure coincides 
with the epoch of model 1490 shown in panel (c) of Fig. [T] In 
the "non-suppressing" model sequence, the ThMPs are superim- 
posed on the ensuing broad Ltot trough; the luminosity of the ter- 
minal model E (whose oscillatory secular modes were all at least 
weakly stable) lied slightly below the total luminosity of model 
1490. In the pulse-suppressing sequence, on the other hand, the 
linear instability extends to roughly the epoch when the lumi- 
nosity reaches again the level of the initial bump. 

4. Discussion 

In accordance with lHarm & Schwarzschildl d 19721) we found the 
secular eigenspectrum to show an oscillatory secular mode go- 
ing unstable at the epoch of the ThMPs setting in. The secular 
mode Oi remained purely oscillatory through the first two cy- 
cles. Starting with the third cycle, Oi unfolded temporarily into 
a pair of purely monotonic modes, both of which were unstable. 
The monotonic-mode interlude took up roughly 46 % of a cycle 
(at cycle 7) with a slow decrease of this percentage as the ThMP 
episode progressed. After cycle 21, Oi remained again purely 
oscillatory up to the last ThMP cycle (number 30). It was always 
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Fig. 6. The periods of the overstable oscillatory secular mode Oi 
shown in Fig.|4]during the ThMP episode as a function of model 
number The almost horizontally running lines result from the 
unfolding and merging of secular modes. At the corresponding 
phases, the periods go formally to infinity. The thick grey line 
indicates the cycle lengths that were derived directly from the 
Lhb variation of the evolutionary computations. The open circles 
indicate the periods as measured at pulse phases with Lg = 0, i.e. 
at 'equilibrium' epochs. 



around the Lne-maximum phases where the oscillatory Oi mode 
unfolded. Interestingly enough, O3, the other oscillatory secu- 
lar mode with eigenvalue unfoldings, showed these to occur at 
around the phases of minimum Lhe- 

Figure |6] shows the periods computed for the mode branch 
Oi as black continuous line. The lengths of the ThMP cycles 
derived directly from the evolutionary computations are plotted 
as a grey line. The sharp horizontal excursions in the figure are 
caused by the unfoldings of the complex conjugate eigenvalues 
and by the merging of two monotonic secular modes where the 
period goes formally to infinity. The agreement between linear 
theory and the cycle lengths deduced from the nonlinear stellar 
evolution is reasonable. The variation of the "linear" period de- 
rived from the secular stability analysis during a ThMP cycle is 
reminiscent of the variation of the pulsation period computed on 
a nonlinearly pulsating star at various phases of a pulsation cy- 
cle (cf. footnote in Sect |3.2t . The open circles denote the periods 
as measured at phases with Lg = 0, i.e. at formal equilibrium 
phases. The maximum deviations of the cycle lengths derived 
from the evolutionary computations and the periods at equilib- 
rium phases do not exceed 40%. The correspondence of the two 
differently computed cycle lengths in Fig. |6] are considered as a 
good indication that a linear secular analysis is appropriate to at 



least quaUtatively understand the cyclic variability of the ThMP 
phase; this point of view is further supported by the results of the 
secular analyses of the pulse-suppressing model sequence which 
was presented in Sect. 13.31 The periods obtained in the pulse- 
suppressed series agree very well with the cycle lengths mea- 
sured in the evolution computations and they are well bracketed 
by the secular analyses of the thermally pulsing models. 

The rest of the discussion section is structured as follows: 
First, we comment on the local criteria for secular instability and 
the difficulty of their use in case of oscillatory secular modes. A 
short digression on the attempts to understand what discrimi- 
nates between monotonic and oscillatory secular modes follows. 
After that, we delve into the properties of the linear eigenmode 
spectrum as we computed it. We close the discussion section 
with checking the quality of the Gabriel-expansion during the 
short-term evolution of the ThMP phase. 

4. 1 . Local stability criteria 

The paper of Schwarzschild & HarmI (Il965h contains a simple 
conceptual model of a nuclear-burning shell in a star The for- 
mulae they obtained at the end of their local analysis allow to 
quickly estimate the shell's stability against thermal perturba- 
tions; therefore, the local criteria were used ever since in com- 
putational research and in textbooks. It turns out that two condi- 
tions must necessarily be fulfilled for instability: The shell must 
be thinner than 



r 4 ■' 



AT 4 

but sufficiently thick so that — > — , (1) 

T et 



the quantity / measures the deviation from a homologous pertur- 
bation; most frequently / is chosen ad hoc to be about unity, i.e. 
close to homology. The variation of physical quantities across 
the nuclear burning shell are indicated with a leading A. The 
width of the shell is not uniquely defined and gives rise to consid- 
erable ambiguity. The remaining quantities take on appropriate 
mean values that are representative of the nuclear-burning shell. 
With St we denote Sloge/cJlog T at constant pressure. The rest 
of the quantities have their canonical meanings. Note that the 
instability ffiat is described with conditions ([TJ is a monotonic 
one; this simple one-zone - type model does not admit of com- 
plex solutions at all. Therefore, relying on criteria ([T]l in the case 
of thermal-pulse cycles that start via oscillatory secular modes 
can be misleading. 

Table[T]contains numerical values of the two sides of both in- 
stability criteria in ([1]) for selected stellar models and for various 
choices of the width of the He-burning shell. The thickness of the 
shell was measured as the extension over which nuclear energy 
generation exceeds a prescribed fraction of Smax- The range of 
numbers given in column is caused by the non-negligible 

variation of sr across the nuclear burning shell. Model 1470 was 
secularly stable, it showed no unstable secular modes and around 
that epoch Lee varied still only weakly and monotonously. At 
model number 1516, Lhb just passed through the first peak and 
mode Oi was overstable at this epoch; the ThMP-phase just de- 
veloped. Finally, model number 3234 lies after the ThMP phase 
with no overstable oscillatory secular modes anymore and Lee 
changed again slowly and monotonously. 

Most frequently, only the first condition of the instability cri- 
teria ([TJ is referred to in the literature; i.e. only the test if a shell 
is sufficiently thin for the instability to develop. The instability 
condition in expressions ([T]i is equivalent to the physical notion 
of a generalized specific heat (also referred to as "gravothermal" 
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Table 1. Evaluation of physical quantities across the He-shell 



Model 


width 


Ar/r 


* 1 / ^ 


Ar/r 


4/et 


1470 


/ 10 

"-■max / ^ ^ 


0.122 


0.41 


0.067 


0.1 - 0.09 




Emax/lOO 


0.23 


0.41 


0.129 


0.1-0.09 




fimax/1000 


0.318 


0.41 


0.168 


0.1-0.08 


1516 


emax/10 


0.175 


0.41 


0.067 


0.11-0.09 




£max/100 


0.258 


0.41 


0.109 


0.11-0.09 




£max/1000 


0.349 


0.41 


0.149 


0.11-0.08 


3234 


emax/10 


0.413 


0.41 


0.149 


0.32-0.11 




Smax/lOO 


0.514 


0.41 


0.206 


0.32-0.13 




£max/1000 


0.585 


0.41 


0.224 


0.32 - 0.08 



specific heat) Q which goes positive within the nuclear burning 
shell. lOiannone & Weigert ( 1967) generalized the local instabil- 
ity criteria ([T]l, adopting a more general equation of state and al- 
lowing for more variation of the physical bac kground quantities 
across the nuclear-burning shell. lYoon et all P004) elaborated 
on the local stability criteria of nuclea r burning shells; building 
upon the approach of lGiannone & Weigert (1967)Flbut still with 
the simplifying assumption of homologous perturbations. The 
aim was to quantitatively test stellar evolution models towards 
thermal-pulse instabilities. The problem with the local stability 
criteria lies mainly in the necessity to have to assume a func- 
tional form of the perturbations (usually homologous motion) 
and therefore to constrain the nature of the instability and to have 
to ado pt a not well definable thickness of the shell. lYoon et akl 
(12004 mention that the shell thickness in their helium-accreting 
white dwarf changes by about 15% if they use once the radii at 
r(e - 10"^£,nax) and once at r{e = lO^-'smax) with s^m being 
the peak nuclear burning rate. For their discussions, Yoon et al. 
chose r(s = 2 x 10"^ ■ Smax) as then their local analyses agreed 
best with the evolutionary computations. This approach might be 
incomplete, however, since, as we have discussed before, the dis- 
sipation inherent in the computational schemes for stellar evolu- 
tion might suppress thermal pulses despite their being physically 
present and detectable in a full linear secular analysis. 

In TablelU local secular instability obtains if the numbers of 
column "Ar/r" are smaller than those of column "ri/4" and if 
the numbers of column "Ar/T" are bigger than those of column 
"4/st"- For model 1470, the geometrical criterion in ([T]) indi- 
cates instability for all choices of shell thickness; the thermal 
part of the local condition favors instability for sufficiently thick 
shells (i.e. for Smax/lOO and Smax/lOOO). Note that for model 
1470 (cf. Fig.lU neither the solution of the full secular boundary- 
value problem nor stellar evolution shows a sign of secular in- 
stability yet. For model 1516, the geometric and the thermal 
local criteria behave the same as for model 1470. In this case, 
it is compatible with stellar evolution and full secular-analysis 



'* To explain the physical processes during a secular instability, 
iKippenhahnI jl970h was the first to refer to an "effective specific heat" 
that goes negative. The baptizing of this effecti ve specific heat as 
"grav othermal specific heat" seems to go back to ISugimoto & Mivaiil 
( 119811) who were guided by an analogy with the gravothermal catastro- 
phe of self-gravitating systems. 

' The^ Giannone & Weigerj i 19671) derivation of the local stability cri- 
teria of nuclear burning shells is more general than what Schwarzschild 
& Harm (1965) (SH65) did. Not only do Giannone & Weigert accom- 
modate a more general equation of state than SH65, they are also not 
resorting a priori to homologous perturbations. The resulting stability 
criteria look considerably more crowded; on the other hand one gains 
a broader range of applicability than the SH65 one. This is a point of 
interest when objects other than low-mass AGB stars are studied. 



results. For model 3234, at an epoch in the post-ThMP phase, 
all choices of shell size point to stability and, depending on the 
choice of a representative magnitude of 4-/st, the same applies 
to the thermal criterion. 

In agreement with Yoon et al.l (|2004|) we find that the par- 
ticular choice of shell thickness influences crucially the result 
of the local criterion. Moreover, the onset of the ThMPs is not 
captured reliably with the local stability criteria. This does not 
further surprise as the ThMPs are the result of the instabiUty of a 
complex secular mode and the local criteria are based on a one- 
zone model description that do not support complex solutions at 
all. 

In the thermal criterion of conditions ([TJ, thickness is ex- 
pressed as temperature and luminosity differences across the 
burning shell, these quantities are more sensitive to a particular 
choice of shell-boundaries than the radius measures; therefore 
the thermal criterion is more sensitive to a particular choice of 
shell thickness. 

4.2. Oscillatory versus monotonic secular solutions 

Aheady iGabriel & Noelsl d 19721) stressed that due to the non- 
hermiticity of the secular stability problem, complex eigenval- 
ues must not be considered unusual. An unambiguous prescrip- 
tion of conditions that must prevail in a star to lead to complex 
eigensolutions have not been put forth in the literature so far 

As mentioned at the beginning of Sect. |4] the merg- 
ing/unfolding phases of the Oi mode branch during the pulse cy- 
cles differ from that of O3 (see Fig.|4]i; this immediately enforces 
the important conclusion that not stellar background quanti- 
ties alone discriminate between monotonic and com plex secula r 
eigensolutions, an aspect already conjectured by Defouwl ( ll973h . 
Despite adopting the point of view that the ther mal pulse s them- 
selves are (nonlinear) relaxation oscillations, iDefouwl (1 19731) 
emphasized the combination of hydrostatic readjustment to a 
heat perturbation and the magnitude of the temperature depen- 
dence of nuclear burning, st, to decide ove r the nature of the 
secular eigenvalue. Earlier on. lDennisl(ll97lh studied the depen- 
dence on St of the secular stability of the helium-burning shell 
of massive stars. Low values of st, left the considered secular 
mode oscillatory damped; increasing st drove it overstable and 
for St > 53, it eventually unfolded into two monotonic branches. 
In our computations, we found the magnitude of st at the peak 
of the unstable He-burning shell to stay at about 40 and to vary 
by less than 5% during the ThMP cycles; in the case of the un- 
folding and re-merging o f the Oi mode, this variation was fur- 
thermore opposite to what' Den"ni si (119711) suggested for this pro- 
cess. Hence, together with the phase-shifted merging/unfolding 
behavior of O3, we conclude that the importance of st lies at 
best in its magnitude but not in its variation. 

4.3. On the secular-eigenmode spectrum 

Looking at the modal diagram in Fig. H] we realize that the os- 
cillatory and the monotonic secular modes evolve differently on 
the (Ti - age plane . This applies also to cyclically repeating 
episodes of unfolding of complex eigenmodes into two mono- 
tonic modes during the ThMP phase and in particular it is true 
for the early ThMP phase when the He shell source takes over 
from central helium burning. However, throughout the whole 
ThMP phase the stable oscillatory secular modes (O2 and O3) 
show an average cri increase that exceeds that of the mono- 
tonic modes; this suggests that oscillatory and monotonic sec- 
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ular modes are each influenced by different regions of the stellar 
interior Unfortunately, the secular eigenproblem admits of no 
special mathematical structure, so that all the helpful tools de- 
veloped for the hermitian adiabatic pulsation problems do not 
grip. We must therefore resort to analyze eigenfunctions at se- 
lected instances and hope to derive therefrom useful hints on the 
nature of the secular modes. 

In the following, we refer to the eigenfunction components ^, 
p, t and 6q = T 6s standing for the relative radius, pressure, tem- 
perature, and heat perturbation, respectively. We studied in more 
detail the eigenfunctions at the begin of the ThMP episode: For 
model 1488, just before the onset of the thermal micropulses, for 
model 1492, one stellar model after the onset of the micropulses, 
and for model 1516, at the Lg - Q phase after the first pulse cy- 
cle. Since, as mentioned before, the secular problem admits of 
no favorable mathematical structure, we cannot assume any or- 
dering of the modes (e.g. no correlation of node numbers and 
the magnitude of the eigenfrequency) and indeed this is what we 
encountered. The counting was exemplarily done on eigenfunc- 
tions of the lowest five monotonic eigenmodes, and the first three 
oscillatory ones (Oi, O2, O3). Node counting revealed that the 
mechanical eigenfunction components and the thermal ones be- 
have differently. For all modes considered, the mechanical eigen- 
function components had either 2 or 3 nodes; i.e. the number of 
nodes of the displacement and the relative pressure perturbation 
was independent of the particular mode. Furthermore, the num- 
ber of nodes of a chosen mode was subject to change (by one 
or two nodes) as the star model evolved. For all secular modes, 
monotonic and oscillatory, the radial displacement admits a node 
either within or close to the hydrogen burning shell. 

Figure I2] plots representative eigenfunctions of mode Oi in 
model 1516. The top panel shows the spatial variation of the 
phase of the heat perturbation 6q. The middle panel displays the 
real parts of t, and of the suitably scaled 6q\ the correspond- 
ing imaginary parts of the eigenfunctions are contained in the 
bottom panel. For better orientation, the middle panel is com- 
plemented with a dash-dotted line tracing the (again suitably 
scaled) spatial variation of the star's luminosity. The node of f at 
the H-burning shell is reminiscent of the "mirror principle" en- 
countered in the spatial variation of the radius evolution of many 
shell-burning stars. Only mode Oi has a relatively large displace- 
ment amplitude also interior to the helium burning shell. The 
displacements of modes O2 and O3 become marginal already 
below the hydrogen-burning shell. On the other hand, many of 
the monotonic modes which we studied have non-negligible dis- 
placement amplitudes in the deep interior with these amplitudes 
growing with |cri|. Also the monotonic secular modes admit of 
a node close to the H-burning shell. The outermost node lies in 
the envelope with its position in radius to grow with |cri|, for 
all these monotonic modes, the maximum amplitude is reached 
at the surface. For all modes considered, the perturbations be- 
come homologous below logjc « -3.5; i.e. homologous vari- 
ation prevails only over a tiny fraction of the star's radius and 
in particular it does not obtain at and above the nuclear burn- 
ing shells. Focusing on the perturbations of thermal quantities 
revealed that e.g. 6q of the oscillatory secular modes (Oi to O3) 
is considerably confined to the intershell region of the stellar 
models with some noticeable amplitude at the base of the star's 
envelope, just above the hydrogen shell. In contrast, the heat per- 
turbations of the monotonic modes achieve their largest ampli- 
tudes in the very deep stellar interior with the heat perturbations 
in the nuclear-burning shell region being relatively small. Based 
on the behavior of the thermal perturbations, we can consider the 
oscillatory modes as functionally different from the monotonic 
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Fig. 7. Spatial run, measured in fractional radius x, of eigenfunc- 
tion data of the Oi mode in model 1516. The two lower panels 
show the real and imaginary parts of the relative temperature per- 
turbation t, the radial displacement ^, and the heat perturbation 
Sq. The heat perturbation and the luminosity L (dash-dotted line 
in the middle panel) of the stellar model are suitably scaled to fit 
the plot range. The top panel shows the spatial variation of the 
complex phase angle of the heat perturbation 6q. 



ones. Furthermore, the thermal eigenfunctions show a node be- 
havior that differs from that of the mechanical ones. The number 
of nodes of the thermal components exceed that of the mechan- 
ical ones in all models considered. In contrast to the mechani- 
cal eigenfunctions, the thermal eigenfunctions of different mode 
branches have always different numbers of nodes (e.g. at model 
1488: 4 nodes in 5q for Oi, 7 for O2, and 1 1 for O3). But again, 
as for the mechanical eigenfunctions, the number of nodes of the 
thermal eigenfunctions of a particular mode branch is not con- 
stant as the model evolved. 

It i s worthwhile to point o u t that , in contrast to the find- 
ings of iHarm & Schwarzschildl (11972 ). the imaginary parts of 
the eigenfunctions of the oscillatory modes are always compara- 
ble to the real parts over most of the stellar interior; this applies 
for the mechanical quantities as well as for the thermal variables 
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(cf. Figure |7]). The phase of 6q varies significantly over most of 
the star's radius, in particular the envelope is not varying syn- 
chronously; a fact which is quantified in the top panel of Fig. |7] 
The most rapid phase variations occur in the H-burning and just 
below the He-burning shell. Across the helium- burning shell, 
the phase stays essentially constant, so that eigenfunctions be- 
have, at least locally, like standing waves. And as to be expected, 
in the region where the perturbations are homologous, i.e. for 
log X < -3.5, the phase variation levels off. 

The phase variation of 6q i s also relevant for another aspect: 
iHarm & Schwarzschildl d 19721) postulated without further speci- 
fication that the complexification of secular eigenfrequencies re- 
lies on the coupling of a "reacting" and a "driving" layer, with 
6q being 7r/2 out of phase in the two layers. The reacting layer 
was identified with a region just inside the He-burning shell; 
stellar material above the helium shell was thought to have too 
low a heat capacity to constitute a pertinent reactive layer and 
w as therefore dis carded. But only the two-zone models (TZM) 
of lDefouwl (Il973h eventually characterized the physical require- 
ments for the occurrence of complex eigenfrequencies in the sec- 
ular problem. The TZM were successful to reproduce results of 
previous numerical experiments that aimed at understanding the 
occurrence of oscillatory secular modes. 

To quantify the phase difference between reacting and driv- 
ing region we resort here to the "ultrasimplified model" which 
negle cts the thermal coupling between the two zones (cf. Sect. 
VI of iDefouwl 1 1 973h ; in Defouw's notation the phase angle (p 
between the heat perturbations of the TZM can be written as 
tan0 = -nil {A\Ci\ — hr). Hence, a phase difference of 7r/2 (as 
referred to in lHarm & Schw arzschild 1972) between the driving 
and the reacting zone is recovered only for the very special cir- 
cumstance of A\Ci\ - riR - 0; i.e. for hr - Q and cn = which 
means for a neutral oscillation (temporal variation oc exp(« ■ t)) of 
a model with an infinite gravothermal specific heat of the driving 
region, or for the very special correlation AiC\\ - ««. In all other 
cases, however, < < njl between 6q of the driving and re- 
acting regions, depending on URlni and A\C\\lni. Together with 
the fact that the heat capacity of the intershell region not being 
lower than that of the He-shell, which is considered as the driv- 
ing region, we conclude that it is not unreasonable to contem- 
plate the intershell region or part thereof as the reacting layer 
Furthermore, the differential work plot shown in Fig.|8]hints at a 
higher reactivity of the intershell region than of the stellar layers 
on the inside of the He shell. 

The lower two panels of Fig. |7]display the heat perturbation 
6q as a heavy line and the relative temperature perturbation, t, 
as a thin solid line. The displayed spatial run is representative 
of the Oi behavior during the whole ThMP phase. We empha- 
size the positive correlation between 6q and t in the He-burning 
shell: Locally a necessary condition for a secular instability to 
develop; if a star goes eventually unstable depends of course on 
the detailed damping and excitation efficiencies throughout the 
whole object. (The amplitude of 5q, an absolute perturbation, 
was divided ad hoc by 2 • 10'^ to match the scale of the fig- 
ure.) \f 6q - A ■ 6T/T with A > lies in nuclear-burning layers 
with sufficiently temperature-sensitive reactions, this helps for 
a marginal heat perturbation to develop into an instability. The 
6q - A ■ 6T /T correlations in suitable stellar layers was encoun- 
tered only for Oi modes. The other oscillatory secular modes, 
despite having also positive Sq - 6T/T correlations, these corre- 
lations were confined to the intershell region where no physical 
agent exists to exploit the temperature perturbation at appropri- 
ate phase. Also the monotonic modes show regions of suitable 
6q - 6T/T correlations which were usually located in the very 




-4 -3 -2 -1 logx 



Fig. 8. Differential work integrals -dw/dlogr throughout a 
star's interior The sign convention is so as to show driving re- 
gions positive and damping ones negative. The dash-dotted line 
depicts, to help orientation, the spatial variation of the total rel- 
ative luminosity, L/L,. The data is shown for model 1492, the 
epoch just after the onset of the thermal micropulses with mode 

01 having gone overstable. 

deep interior, below the He-burning shell, i.e. in the inert C/O 
interior 

To more directly illustrate the action of the helium shell to 
power the instability and thermal activity of the intershell re- 
gion we computed work integrals for some overstable oscilla- 
to ry secular mode s. Following the reasoning and interpretations 
of iGlatzell d 19941) . we justify the adoption of the work integral 
even in the case of rapidly growing instabilities. 

For model 1492, Fig. [8] displays the differential work curves 

dw 1 

-— oclm(t*p)Pr^ 6, (2) 

dlogr 

of the oscillatory modes Oi to O3. We defined 6 = 
d\ogpld\ogP at constant temperature. The asterisk refers to 
the complex conjugate of the quantity to which it is applied. All 
differential-work curves in Fig. |8] are normalized independently 
to unity at the positions of their respective maximum driving as 
we are mainly interested to illustrate the domains of "activity". 
The spatial variation of the differential work is representative of 
each mode for the whole ThMP episode; it is only the strengths 
of the various driving and damping regions that change relative 
to each other For better orientation, the thin dash-dotted line in 
Fig. |8] traces out the spatial run of the total relative luminosity. 
For the Oi mode, it is evident that the maximum driving oc- 
curs in the helium-burning shell. The intershell domain between 
helium and hydrogen burning constitutes the dominant damp- 
ing region. Some minor additional damping is observed in the 
outermost layers. As the ThMPs begin, the driving in the He- 
burning shell overcompensates the intershell damping. Finally, 
as the star stabilizes secularly, the driving and damping features 
persist as before; however, the damping regions eventually win 
over the weakening driving effect in the fattening He shell. For 

02 and O3 , only weak driving occurs in the He-burning shell, 
the intershell region exhibits almost canceling driving and damp- 
ing. It is then the envelope that eventually determines the fate of 
the modes' stability with strong damping throughout most of the 
out er layers . 

iHansenl ( Il972h pointed out the running thermal-wave char- 
acter of oscillatory secular modes. Studying the temporal varia- 



A.Gautschy and L. G. Althaus: Oscillatory Secular Modes 



11 



Ar/r 



0.005 - 



0.000 _ 



-0.005 - 



-0.010 - 




-4 -3 -2 -1 

logx 

Fig. 9. The top panel shows the relative radius change, Ar/r, as 
a function of the relative stellar radius x derived from the stellar 
evolutionary sequence, as the star evolves through the ThMP- 
cycle between model numbers 1732 and 1814. This phase covers 
the largest-amplitude mode-unfolding episode depicted in Fig.|4] 
The bottom diagram compares the quantity Ar/r at model 1792 
with the result of the Gabriel-expansion with 9 terms, applied to 
the displacement eigenfunction. 



tion of the complex eigenmodes in our models revealed indeed 
a pronounced propagative behavior. All complex eigensolutions 
behave similarly in that outward propagating patterns develop in 
the inter-shell region; the He-burning shell acts as a - spatially 
stationary - piston for these features. Below the He shell, heat 
perturbation and radial displacement behave mostly like stand- 
ing waves with only a weak inward-running component. 

4.4. The relevance of the secular modes on the star's 
evolution 

At the end, we want to find out how well the iGabriell dl972h - 
expansion does in the context of our ThMP problem. We ad- 
dress this aspect with the help of Fig. |9j its top panel shows the 
spatial variation of Ar/r = 1 - r{m,t())/r{m,t\) (with fo, t\ de- 
noting temporally neighboring models) during the ThMP cycle 
which is bounded by models 1732 and 1814-this cycle encom- 
passes the largest-amplitude mode-unfolding seen in Fig.m The 
model numbers given in Fig. |9]refer to the later epochs used to 
compute Ar/r; the bump at about logx = -2.8 (location of the 
He-shell) is spatially an essentially stationary feature, the same 
applies to the homologous radius variation of the even deeper ly- 
ing regions. An almost perfect node (at logx -2.2) lies at the 
base of the hydrogen-burning shell. This node is a manifestation 
of the "mirror principle" which is referred to frequently in the 
context of shell burning stars. Finally, an outward propagating 
bump develops above the hydrogen-burning shell; its propaga- 
tion into the star's envelope can be traced even in the time-series 
of Ar/r in the top panel of Fig.|9] 



IGabriell ( Il972h showed that for sufficiently short time-scales, 
over which a star's evolutionary variations can be linearly ap- 
proximated relative to a reference epoch, the physical quanti- 
ties can be expanded into a series of secular eigenfunctions (see 
eq. (13) in IGabriell 19721) . We followed Gabriel's procedure to 
find out how well the linear eigensolutions approximate the vari- 
ation induced by the star's short-term evolution. At (between 100 
and 300) equally spaced gridpoints in log x, we solved the least- 
square problem for the coefficien ts of the series expansion in 
T5s (cf. eq. (15) in IGabriell [T9T2 ). Table |2| lists the results for 
model 1792 for various choices of terms included in the series 
(first column). For each number of modes included in the series 
expansion, the sum of the coefficients was normalized to unity 
separately. The mode designated by O3 was unfolded at model 
1792 so that its expansion coefficient is real and not complex 
as for Oi and O2. The mono tonic modes which were included 
in the series expansion are labeled as M,. In contrast to Gabriel's 
gravitationally contracting model stars, our ThMP models do not 
show an equally rapid convergence of the series as a function of 
the number of modes (i.e. terms) included. In accordance with 
the eigenmode properties discussed in Sect. l4.3l we found modes 
Oi and O2 to be essential to approximate well the functional be- 
havior in the nuclear burning shells; to also match the core and 
the envelope variation, the inclusion in the series expansion of 
(damped) monotonic modes, in particular of Mi, turned out to 
be indispensable; however, at the price of added spatial small- 
amplitude oscillations in -4 < logx < -3. An impression of the 
quality of the result of Gabriel's procedure applied to Ar/r of 
model 1792 with 9 terms in the series expansion can be got from 
the bottom panel of Fig.|9l 

The Gabriel ( 1972) - expansion turned out to do reasonably 
well to approximate the spatial variation of perturbed physical 
quantities (see bottom panel of Fig. O. However, in contrast to 
many stellar pulsation problems, linear secular theory failed to 
reproduce the amplitude evolution of the variability throughout 
the ThMP episode. This is not further surprising as this aim 
lies outside the domain of validity of the Gabriel-expansion. 
Nevertheless, the period of the overstable secular mode during 
the ThMP phase provides a good approximation to the pulse- 
cycle length of the fully nonlinear problem. 

Finally, it is worth emphasizing that the secular problem 
is strongly non-normal. So that we cannot, in general, expect 
a close match of the linear results with the nonlinear develop- 
ments. For non-normal operators, considerable transient phases 
are usually encountered at the onset of instabilities, so that 
normal-mo de eigenanalyses might have only limited pr edictive 
power (e.g iTrefethen & Embreg|2005l ; lEisenma nl l2005l) . In our 
case here, though, the agreement is surprisingly good. 

5. Conclusions 

A large number of 3M0 stellar models that exhibit thermal mi- 
cropulses at the onset of their double-shell nuclear burning phase 
were studied by means of linear secular eigenanalyses. The 
micropulses are considered as physically comparable with the 
larger-amplitude thermal pulses of low- and intermediate-mass 
stars usually encountered along the advanced AGB. Our compu- 
tations showed that the ThMP episode sets in via the instability 
and finally dies away via the restabilization of an oscillatory sec- 
ular eigenmode. And most importantly, the overstability of the 
lowest-frequency oscillatory secular mode persists over most of 
the ThMP phase. During the numerous ThMP cycles, the oscilla- 
tory modes unfold temporarily into pairs of purely monotonous 
secular modes. This phenomenon, we think, is the result of per- 



12 



A.Gautschy and L. G. Althaus: Oscillatory Secular Modes 



Table 2. Coefficients of the expansion of rd,s|e,.o/ in terms of secular eigenmodes for model 1792. The first column gives number 
of terms of the series expansion. The modes used in the series are denoted by O's for the oscillatory ones (cf. Fig. |4]i and M's 
for monotonic secular modes. The second lines in the columns of Oi and O2 give the imaginary parts of the respective expansion 
coefficients. 



number of modes 


Oi 


O2 


O3 


M, 


M2 


M3 


M4 


M5 




3 


0.177 
0.427 


0.301 
0.124 


-0.031 














4 


0.238 
0.513 


0.355 
0.270 


-0.085 


-0.289 












5 


0.223 
0.494 


0.331 
0.301 


-0.074 


-0.289 


0.003 










6 


0.264 
0.536 


0.363 
0.236 


-0.087 


-0.319 


0.012 


-0.006 








7 


0.248 
0.456 


0.331 
0.360 


-0.049 


-0.293 


0.014 


-0.012 


-0.034 






8 


0.209 
0.387 


0.105 
0.474 


-0.012 


-0.278 


0.012 


-0.009 


-0.012 


0.126 




9 


0.208 
0.383 


0.158 
0.481 


-0.009 


-0.274 


0.012 


-0.009 


-0.016 


0.084 


-0.016 



forming stability analyses on the actually pulsing stellar mod- 
els. We conclude that the ThMPs are the nonlinear development 
of the linear instability of the oscillatory 0\ mode. Most im- 
portantly, the cyclic nature of the ThMPs has its origin in an 
unstable linear oscillatory secular mode that persists throughout 
the ThMP episode; it is therefore not surprising that the linear 
periods (cf. Fig. |6] and in particular Fig. |5]l correlate w ell with 
the l engths of the nonlinear pulse cycles. The view of Defouwl 
( 1 1973.) of the thermal pulses to be nonlinear relaxation oscilla- 
tions of monotonic type appears to have been premature as it re- 
lied probably on the insufficiently long simulation computations 
of[H arm & Schwarzschild ( 1972). Would the latter authors have 
had at their disposal data covering more thermal-pulse cycles, 
we are confident that they would have observed a re-merging of 
the two monotonic mode branches into a complex pair during 
later phases of pulse cycles with sufficiently low amplitude (cf. 
Fig. |4ll. The linear aspect of the instability which was empha- 
sized in this paper is not capable to capture the nonlinear, finite- 
amplitude development of the various physical quantities during 
the ThMPs. Indications that nonlinear processes are indeed im- 
portant are hinted at in Fig. [3] which shows that the amplitudes 
of the radius and density variations in the nuclear-burning shell 
centers saturate after about the first four ThMP cycles and in 
Fig. |4] with the amplitude of Lhb decreasing slowly after about 
397 My despite the persistent overstability of mode Oi until 
about 403 My. 

Applying the explanation of thermal pulses as the nonlin- 
ear development of a linearly overstable secular mode leads 
also to a natural explanation of the origin of the core-mass - 
interfl ash-period relatio n which was derived along the advanced 
AGB dPacz vnskil 1 975l) . We conjecture that the interflash-period 
is the period of the overstable secular mode which gets shorter 
with increasing core mass. Unfortunately, the evolutionary phase 
of the ThMPs is so short that the core mass does not change 
significantly; therefore, our conjecture could not yet be under- 
pinned quantitatively. Nevertheless, Fig. 9 hints at a decreasing 
period of the overstable secular mode as the star model ascends 
the AGB. A detailed study of thermally pulsing AGB stars re- 
mains to be done. 

Well known, but frequently not paid sufficient attention to is 
the fact that stellar-evolution results depend, partly even sensi- 
tively, on the time-stepping of their computation. As known for 



a long time (e.gl Rosdl 19661) . we could essentially suppress the 
pulses when adopting even a moderately large time-step (a factor 
less than two compared with the pulse-resolving time-step). The 
nonlinear "pulse-suppressing" evolution computation showed a 
very weak signature in Lhc of three strongly damped thermal 
pulses only. The linear stability analysis on these models contin- 
ued, on the other hand, to reveal an unstable oscillatory secular 
mode with the instability extending over essentially the same 
evolutionary phase as found in the "non-suppressing" computa- 
tion. Hence, the stellar evolution equations as discretized in most 
computer codes are inherently dissipative and they are prone to 
suppressing secular instabilities at large time-steps. Therefore, 
depending on the particular goal, it is not advisable to eagerly 
aim at evolving stars with a minimum of computational time- 
steps; it also seems dangerous to estimate the stability of shell 
sources via local criteria alone in stellar-evolution computations 
(cf. Sect.gl]!. 

The possibility to decide observationally if ThMPs are real- 
ized in nature seems bleak; the luminosity variations during the 
pulses are very small and on much too long a time-scale to be 
potentially detectable in flagranti. Furthermore, the number of 
stars in the appropriate evolutionary stage in a stellar aggregate 
such as a cluster seems too small to be observable via an accu- 
mulation effect in say a color-magnitude diagram. Last but not 
least, we found no evidence of interpulse mixing in the envelopes 
of ThMP-ing stars so that spectroscopic forensics to identify 
such a phase appears implausible. Nonetheless, from the point 
of view of understanding stars, the ThMPs constitute a benefi- 
cial computational laboratory to study fundamental aspects of 
secular stability of complex stellar models in advanced evolu- 
tionary stages pertinent also to the nature of thermally pulsing 
upper- AGB stars. 
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Appendix A: Approximate stability analysis of stars 
in thermal imbalance 

In the following, the linear equations are compiled as they were 
used to compute the radial secular modes of stars in thermal im- 
balance. We relied on the perturbed Lagrangian quantities 



6r 6P 6T , 6L 

^^7' P^T^ 



(A.l) 



The time dependence of the problem is factorized by exp(icrf). 
The eigenvalue of the resulting linear problem, cr, is normalized 

by \\cr\\ - ^J3GMt/Rl. The linearized stellar structure equations 
for radiative regions, including approximated thermal-imbalance 
terms, are cast into the form 
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with the secondary quantities in the matrix defined as 

523 = Ci (C2 Sp + Tp+ icrVad) , 

siA - Ci (C2 bt + Tt - icr) , 

S3i = (4- + 3ciCr^)V , 
S44 = V Vi-ad (4 - Xt) , 

with ci - (r/Rt)^ /(Mr/Mt) and x - r/Rt. Furthermore, 
MrUcpT\\o-\\ 



Ci = 
C2 = 
C3 = 



cpT\\cr\\ 
d, log T 



C4 = V, 



Ml 
d, logP 
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The quantities U and V represent the two homology invariants. 
The temporal derivatives of log P and log T are the ones used in 
the gravitational-energy term of the stellar-evolution computa- 
tions. Finally, we approximate the thermal-imbalance terms as 



Tp - -Cjcpp + C4i6p — a) 
Tj — -CjCpT + C4 (St + S) . 



(A. 14) 
(A. 15) 



With cpq we refer to dlogcp/dlogq, the same applies to 
6q. The remaining not explicitly defined variables have ei- 
ther their ca nonical meaning as used in stellar structure the- 
ory (see e.g. Kippenhahn & Weiger3 1 19941) or as referred to in 
iGautschv & Glatzell diggoT Regularity constraints in the star's 
center enforce the central boundary conditions; at th e stellar sur- 
face, t he physically motivated ones used in iGautschv & Glatzell 
([T990h were implemented. 
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